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Anderson localization has been observed for a variety of media, including ultracold atomic gases 
with speckle disorder in one and three dimensions. However, observation of Anderson localization 
in a two-dimensional geometry for ultracold gases has been elusive. We show that a cause of 
this difficulty is the relatively high percolation threshold of a speckle potential in two dimensions, 
resulting in strong classical localization. We propose a realistic point-like disorder potential that 
circumvents this percolation limit with localization lengths that are experimentally observable. The 
percolation threshold is evaluated for experimentally realistic parameters, and a regime of negligible 
classical trapping is identified. Localization lengths are determined via scaling theory, using both 
exact scattering cross sections and the Born approximation, and by direct simulation of the time- 
dependent Schrodinger equation. We show that the Born approximation can underestimate the 
localization length by four orders of magnitude at low energies, while exact cross sections and scaling 
theory provide an upper bound. Achievable experimental parameters for observing localization in 
this system are proposed. 


I. INTRODUCTION 

Anderson localization (AL) [1, 2]—the preclusion 
of wave propagation in a disordered medium by 
interference—has been observed in many settings, 
including, e.g., for light and sound [3, 4]. Recent ex¬ 
periments observing AL for ultracold atomic gases 
expanding in disordered optical speckle potentials 
show promise for gaining enhanced understanding 
of how microscopic disorder characteristics affect lo¬ 
calization and the interplay with inter-particle in¬ 
teractions. In particular, these systems allow for in¬ 
dependent control of inter-atomic interactions and 
the disorder strength and correlation length. Lo¬ 
calization has been observed and its dependance on 
these disorder parameters studied for gases confined 
to one-dimensional geometries [5, 6] and in three di¬ 
mensions [7-9]. In two-dimensional gases, however, 
the classical diffusive regime [10] and the impact of 
disorder on superfluids have been explored [11, 12], 
but AL has not yet been observed. Studying lo¬ 
calization in two dimensions using ultracold gases 
is especially desirable given the many outstanding 
questions regarding localization in two-dimensional 
electronic solids [13-15]. 

In this paper, we discuss how classical trapping 
effects complicate the observation of AL for speckle 
disorder in two dimensions. We demonstrate how 
point-like disorder [Fig. 1] avoids these problems. 
Through a combination of analytical and numeri¬ 
cal simulation we show that a disordered potential 
of this type would enable experimental observation 
of two-dimensional AL for ultracold gases. Point¬ 
like, two-dimensional disorder and time-dependent 
simulations of localization have not been considered 
in previous theoretical studies of ultracold systems 


[16-18]. Like previous studies, we consider an ideal 
two-dimensional geometry and ignore inter-particle 
interactions and quantum statistics. As our focus is 
on realistic experimental conditions, we have limited 
our investigation to experimentally accessible disor¬ 
der strengths and energy scales for ultracold, spin- 
polarized atoms. We expect that this technique 
can be straightforwardly extended to other experi¬ 
mental configurations. 

This paper is organized as follows: in Sec. II 
we describe classical trapping in two-dimensional 
speckle potentials using percolation theory, and we 
show how point-like disorder can avoid this problem 
for realistic experimental parameters. In Sec. Ill we 
identify experimentally accessible parameters for lo¬ 
calization in point-like disorder using scaling theory 
and the Born approximation, which we show may 
fail in this regime. Therefore, in Sec. IV, we calcu¬ 
late exact differential cross sections using numerical 
simulations to more accurately determine localiza¬ 
tion lengths via scaling theory. We also use exact 
time-dependent simulations of wavepacket propaga¬ 
tion in point-like disorder to determine localization 
lengths as would be observed in an experiment that 
allows a gas to expand into a disordered potential. 


II. PERCOLATION IN SPECKLE AND 
POINT-LIKE DISORDER 

In this section, we examine the classical trapping 
properties of speckle and point-like disorder poten¬ 
tials in two dimensions using percolation theory, 
which is the study of the random growth of inter¬ 
connected regions in networks. These regions grow 
in size as a parameter is varied until a transition at 
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FIG. 1. Schematic representation of experimental im¬ 
plementation of point-like disorder. Ultracold atoms 
are confined to a quasi two-dimensional geometry us¬ 
ing a sheet of far-detuned light (not shown). The disor¬ 
dered potential is generated by an additional laser beam 
(green) that passes through a holographic optic (disc) 
and is focused on the atoms (red spheres). The atoms 
experience a disorder potential consisting of a random 
arrangement of Gaussian barriers. 


the percolation threshold, at which a connected re¬ 
gion spanning the system is formed [19, 20]. In the 
context of atom transport in a two-dimensional dis¬ 
ordered potential, the problem is easily visualized 
in this way: the disorder potential forms a poten¬ 
tial landscape that constrains the atomic motion. 
Atoms with kinetic energy higher than the perco¬ 
lation threshold will travel freely through the po¬ 
tential, but below a critical fraction of the average 
potential energy atoms will be trapped in potential 
minima of finite size [21]. If this percolation thresh¬ 
old occurs near the same energy scale as Anderson 
localization, AL may not be detectable or experi¬ 
mentally distinguishable from classical trapping. 

Because of their ease of creation and simple sta¬ 
tistical properties, optical speckle fields [22] have 
been used to generate disorder in virtually all ultra¬ 
cold atom experiments exploring disorder-induced 
effects, see, e.g.. Refs. 5, 7, 8, 10-12, 23, 24. A 
notable exception is disorder generated by impurity 
atoms [25]. Optical speckle is produced by focusing a 
laser beam that has passed through a diffuser. The 
atoms experience a potential energy shift propor¬ 


tional to the optical intensity of this light. Speckle 
disorder has an exponential distribution of poten¬ 
tial energies and a spatial autocorrelation that is 
approximately Gaussian. While in three dimensions 
the percolation threshold for a blue-detuned speckle 
field is negligible [26, 27], potential minima appear 
in one and two dimensions that can classically trap 
the atoms. As a result, the percolation characteris¬ 
tics of a speckle potential become a significant con¬ 
straint on attempts to observe Anderson localization 
in reduced dimensions. 

This limitation is in part due to scaling symme¬ 
tries unique to speckle disorder. A distinctive fea¬ 
ture of a speckle field it that it has only one ad¬ 
justable length scale, the correlation length This 
characteristic leads to a simple scaling symmetry 
of the field: any change in scale of the system is 
equivalent to a suitable change in The percola¬ 
tion threshold is necessarily scale invariant, and as 
a result it must be independent of Thus, the crit¬ 
ical energy for a percolation transition in a speckle 
field is always at a fixed fraction of the average 
potential energy A, which simulation shows to be 
roughly 52% in two dimensions [27-29]. This de¬ 
pendence complicates observing AL in two dimen¬ 
sions using ultracold atoms. In infinitely sized two- 
dimensional system, infinitesimal disorder will local¬ 
ize atoms with any kinetic energy [30]. The local¬ 
ization length grows exponentially with the particle 
energy [31], and thus relatively strong disorder is 
required to localize atoms on experimentally acces¬ 
sible length scales. This leads to classical trapping 
of a wide range of particle energies, which may be 
difficult to separate from AL. 

In light of this high, fixed percolation threshold, 
an alternative form of disorder that does not cause 
classical trapping is desirable. The point-like dis¬ 
order we investigate consists of individual Gaussian 
potential barriers with randomly distributed loca¬ 
tions Xi'. 

U(f) (1) 

where Vq is the peak disorder energy, i indexes the 
individual Gaussian potentials, and wj is the rms 
width of an individual barrier. This type of dis¬ 
ordered potential was chosen for its simplicity and 
straightforward experimental realization in a cold 
atom experiment using holographic techniques [32- 
34]. As shown in Fig. 1, a random array of blue- 
detuned focused Gaussian laser beams can produce 
potential barriers as in Eq. 1. The advantage of 
point-like disorder over a speckle potential is that 
the freedom to tune the density n of potential barri¬ 
ers introduces a second length scale and vary¬ 

ing the ratio of w and the average distance between 
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scattering sites allows the percolation thresh¬ 

old to be tuned. 

The impact of this tunability is shown visually in 
Fig. 2. Potential landscapes are displayed for speckle 
and point-like disorder with the same average poten¬ 
tial energy. While the percolation threshold for this 
finite-sized realization of speckle disorder is 0.39A, 
the point-like disorder parameters were chosen to set 
the percolation threshold for the realization shown 
in Fig. 2b to approximately 0.06A. Thus, particles 
with small kinetic energies compared with the dis¬ 
order energy are free to propagate in the point-like 
disorder, while particles with relatively high kinetic 
energies are classically trapped by the speckle disor¬ 
der. 



ily close to zero, and it monotonically increases with 
nw^ towards the limit of the average disorder energy. 
Our interest is in the regime where the percolation 
threshold is lower than that of a speckle potential 
and ideally negligible, while maintaining sufficient 
density of scattering sites such that the particles are 
scattered many times within the system size. 
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FIG. 2. Comparison of percolation in a speckle potential 
(a) and a sparse point-like disordered potential (b). The 
disorder potential energy is shown in false color. The 
colorbar shows the potential energy in nnits of A. Re¬ 
gions in grayscale correspond to energies less than 30% 
the average disorder energy. A classical particle with this 
energy would be trapped in a finite-sized region in the 
speckle potential, but is able to propagate indefinitely 
for the point-like disorder case. 

We used a standard technique to calculate the de- 
pendance of the percolation threshold on the point¬ 
like disorder parameters. Disorder potentials were 
numerically simulated, and the percolation thresh¬ 
old was determined by detecting the formation of 
a connected region spanning the simulation space 
by points below a threshold potential energy. The 
percolation threshold Eth was determined by aver¬ 
aging over independent realizations of the disorder 
potential. While only approximating the percolation 
threshold, which is defined in limit of an infinite sys¬ 
tem, this method gives excellent agreement with the 
known value for speckle disorder [Fig. 3]. In our 
simulation, a fixed system size L is used and the 
number of Gaussian potential barriers and w were 
independently varied. 

As shown in Fig. 3, the percolation threshold was 
found to be a function of only the dimensionless com¬ 
bination nui^, as expected from simple space-hlling 
considerations [35]. For sufficiently dilute disorder 
(i.e., low nw'^), the percolation threshold is arbitrar¬ 


FIG. 3. Percolation threshold for point-like (blue 
squares) and speckle (red circles) disorder potentials. 
The dashed line is the known percolation threshold for 
speckle disorder, and the points are results of our sim¬ 
ulation. For point-like disorder, n and w are indepen¬ 
dently varied, and for speckle disorder the correlation 
length was changed. The inset shows the onset of sig- 
nihcant percolation for point-like disorder. The mini¬ 
mum detectable threshold level is 5 x 10“® in our sim¬ 
ulation. The error bars show the standard error of the 
mean for the average taken over 50 disorder realizations 
at each point-like disorder point and 8 realizations for 
each speckle disorder case. 

We choose to concentrate on nw'^ < 0.03, which 
fulfills both criteria. Under this condition, the per¬ 
colation threshold is less than 2.5 x lO^'^A, which is 
smaller than the percolation threshold for 3D speckle 
[27]. The system size and w are limited by exper¬ 
imental constraints such as optical power, imaging 
signal-to-noise ratio, and numerical aperture. For 
the rest of this paper, we choose to use an experimen¬ 
tally feasible disorder size w = 400 nm. When dis¬ 
cussing proposed experiments, we assume that the 
disorder potential and the light used to create it are 
limited to a 100 x 100 /rm^ area in order to estimate 
the requisite laser power. We use n = 0.2 
(corresponding to nw^ = 0.03) for comparing the lo¬ 
calization lengths predicted by scaling theory using 
either the first Born approximation or the exact scat¬ 
tering differential cross section and for predicting a 
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thermally averaged density profile. In these cases, 
we choose the highest possible n that avoids signifi¬ 
cant classical trapping in order to make AL effects as 
robust as possible. For comparing the results of scal¬ 
ing theory and a simulation of the time-dependent 
Schrodinger equation, we use n = 0.08 /im“^ (cor¬ 
responding to nw"^ = 0.013), which is the largest n 
compatible with our computational resources. For 
n = 0.2 ^m“^ and n = 0.08 /rm“^, there are 
2000 and 800 disorder peaks within the proposed 
100 X 100 area, respectively. 

While we have shown that point-like disorder can 
avoid classical trapping, it must also lead to observ¬ 
able localization lengths to be a viable experimen¬ 
tal option. In principle, atoms in an infinite two- 
dimensional system are localized by infinitesimal dis¬ 
order [31], but only localization lengths smaller than 
the system size are physically meaningful and ob¬ 
servable. Therefore, the rest of this paper is con¬ 
cerned with estimation of the localization length for 
point-like disorder and its dependance on the disor¬ 
der properties and particle energy. 


III. LOCALIZATION LENGTHS: THE 
BORN APPROXIMATION 

In this section we determine analytic expressions 
for the localization length in the limit where the scat¬ 
tered wavefunction amplitude is small. Although 
we will show that this approach is inaccurate in the 
regime of interest, it is useful for developing physi¬ 
cal insight. Analysis of point-like disorder is partic¬ 
ularly simple in this regime, because the differential 
cross section for scattering from a single disorder 
barrier can be determined using the first-order Born 
approximation. This differential cross section can be 
used to determine transport properties in a poten¬ 
tial consisting of many Gaussian potentials, provided 
that there little spatial overlap. Our procedure has 
three steps: we first find the Born approximation 
for the differential scattering cross section da/dO of 
a single potential barrier. We use da/dO to deter¬ 
mine the Boltzmann transport mean free path Ib, 
which characterizes the diffusive properties for par¬ 
ticle transport in the disordered potential. Finally, 
we use scaling theory to estimate the localization 
length ^ from Ib- 

In two dimensions, we write the (unnormalized) 
scattered wavefunction tp{x) as 

= ( 2 ) 

Vr 

where r is the radial coordinate, k is the incoming 
wavevector, / (9) is the scattering amplitude, and 9 
is the scattering angle. We assume a free-particle 


dispersion so that k = ^/2mek/h, where m is the 
mass of the particle, is the kinetic energy of the 
atom, and h = 27rh is Planck’s constant. Given the 
scattering amplitude in the Born approximation 

/ (9) = - (x) d^x (3) 


(f is a unit vector that points in the scattered direc¬ 
tion), the differential cross section for a single Gaus¬ 
sian potential V{x) = Voe~^ is 


da 


\fi0)f 


IT 

Sk 




( 4 ) 


where 9 is the scattering angle. We calculate the 
Boltzmann mean free path Ib, the distance over 
which the direction of momentum remains correlated 
[16], using Eq. 4 for the differential cross section and 
the density of scattering sites n: 


Ib = Is 




( 5 ) 


where the elastic mean free path is G = (na) 
and the total cross section is cr = {da/d9) d9. 
Gombining equations 4 and 5, we find 


Ib = 


1 1 
y/2Tr^ nw 


/ h^k'^ 1 

\^Vo 


2 


( 6 ) 


for the Boltzmann mean free path. 

The connection between Ib and the localization 
length ^ is given in the {klB)~^ ^ 1, L —)• oo limit 
by scaling theory [31] as 

e = (7) 


Thus, given the properties of the disorder and atoms, 
we may calculate a localization length valid in these 
limits using the Born approximation. We choose to 
focus on an experimentally realistic situation. In 
an experiment in which the disorder is generated 
by the dipole force from a far-detuned laser, limited 
laser power constrains the maximum average poten¬ 
tial energy, rather than the value of Vq for individual 
potential barriers. Therefore, we investigate the case 
in which A is fixed, but the density of disorder sites 
n and Vq vary inversely. 

To get a concrete sense of these predictions, for 
disorder generated by a 2 W, 532 nm laser fo¬ 
cused to a Gaussian envelope with a 170 waist 
(which has been employed in experiments on 3D AL 
[7]), realistic disorder parameters are w = 400 nm 
and average disorder strength A = /cs x 1000 nK. 
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The peak potential Vq is determined by the relation 
Vo — Afirnw'^. The resulting localization lengths, as 
a function of particle energy ej, and n, are shown in 
Fig. 4 along with the percolation threshold. For the 
contour lines shown in Fig. 4, the value of Mb varies 
from 0.05-0.3. A localization length of 100 ^m, 
which would result in high imaging signal-to-noise 
ratio and minimal effects from the disorder enve¬ 
lope, can be well-separated from classical trapping 
when there is less than approximately one potential 
barrier per square micron. Furthermore, a regime 
exists for n < 0.25/rm“^ in which classical trapping 
is irrelevant and small localization lengths exist for 
experimentally accessible particle energies. 



FIG. 4. Localization lengths according to scaling theory 
within the Born approximation. Contour lines for three 
localization lengths are displayed using the experimental 
parameters described in the main text. The percolation 
threshold in terms of e*, is shown using solid circles. If 
tk > Eth, then AL will not be observable and the atoms 
will be classically trapped. Because the average disor¬ 
der potential energy is fixed at ks x 1000 nK, Vo varies 
according to Vo = A:b x (l989/n[pm“^]) nK. 


While this standard approach to determining lo¬ 
calization lengths suggests that two-dimensional lo¬ 
calization may be observable using ultracold atoms 
and point-like disorder, it is unclear that the Born 
approximation will be valid in an experiment. Given 
the parameters explored in Fig. 4, accessible local¬ 
ization lengths necessarily involve Vq > ek, and thus 
the Born approximation is suspect. The precise lim¬ 
its of the Born approximation are not always obvious 
[36]. Unlike the three-dimensional case, the Born 


approximation is never valid in two dimensions as 
/c —>■ 0. In two dimensions and for kw ^ 1, a simple 
analysis suggests that the Born approximation may 
be valid for 


Amw Vo 
^ k 


< 1 . 


( 8 ) 


The smallest value for g « 7.5 occurs in the up¬ 
per right hand corner of Fig. 4, and thus the Born 
approximation is not satisfied for this range of pa¬ 
rameters. Similar considerations hold for the high 
energy (i.e., kw ^ 1) regime. Thus, while the Born 
approximation can supply physical insight, the lo¬ 
calization lengths computed in this section may not 
be accurate. 

Using the Born approximation will generally un¬ 
derestimate the localization length. The differential 
cross section in the Born approximation [Eq. (4)] 
depends quadratically on the strength Vq of a scat¬ 
tering site, which accounts for the sharp decrease 
in localization length in Fig. 4 near n = 0, where 
the disorder potential is concentrated in few, very 
strongly scattering sites. However, once Vo is much 
greater than the particle kinetic energy, further in¬ 
creasing Vo must have diminishing effects on scat¬ 
tering. Underestimating the localization length is 
a critical problem because large localization lengths 
may not be observable in an experiment. We deter¬ 
mine the precise deviation from the Born approxi¬ 
mation behavior using numerical simulation. As we 
will show in the next section, the Born approxima¬ 
tion and scaling theory as applied here can under¬ 
estimate the localization length by more than three 
orders of magnitude in the regime we propose to ex¬ 
plore experimentally. 


IV. BEYOND THE BORN 
APPROXIMATION 

In this section we implement numerical simula¬ 
tions of scattering and localization, and use them 
to determine improved estimates of localization 
lengths. We rely on two complementary approaches. 
The first, simplest approach is a single scattering 
simulation: we numerically determine the exact dif¬ 
ferential cross-section for a single scattering event, 
and then use this in place of the analytic result 
(Eq. 4) to compute Ib and ^ as in Sec. III. While 
the values obtained using this method remain valid 
only for Mb ^ 1, they have two important advan¬ 
tages over the analytic solutions: in practice the er¬ 
ror is much less severe at low energy, and, critically, 
the assumptions that entered into the scaling theory 
derivation are such that it provides an upper bound 
on the localization length [31]. 




6 


To apply this method, we numerically extract 
da/dO by a straightforward application of scatter¬ 
ing theory (Fig. 5). We first simulate, via the split- 
step Fourier method [37], a plane wave scattering 
from a single potential barrier (see Appendix A). 
A wavefunction = 0) = jL with energy 

Cfc = /^m and wavevector ky is prepared at ini¬ 

tial time t = 0 in the simulation space, which has 
sides of length L. The wavefunction is propagated 
forward in timesteps 6t according to 

V'(:r,t + dt) = . 

j-i ^ (9) 

where J- represents a Fourier transform. The 
timestep St is chosen to minimize numerical instabil¬ 
ity and such that St ■ V/H <C 1 and St ■ hk"^/2m <C 1. 

While propagating the wavefunction forward 
in time, the Gaussian potential V(x, t) = 
(l — e-‘/o. 2 Tj jg slowly ramped on. We 
determined that the time over which the barrier was 
turned on did not significantly change the simula¬ 
tion results. After a total propagation time from 
r = 1-7 ms (depending on the group velocity), a 
wavefunction distorted by scattering from the po¬ 
tential barrier is produced, as shown in Fig. 5a. For 
the Vq and chosen in Fig. 5, the simulated wave 
is highly distorted, signaling a violation of the Born 
approximation. 

The scattered wave V’sc = '/’{x,t) — tj}'{x,t) is re¬ 
covered using the unscattered wave ip'{x^t), which 
is determined by repeating the simulation with 
V{x,t) = 0. The scattering amplitude fr{d) at ra¬ 
dius r is computing according to Eq. 2 as fr{S) = 
^Jr'i/scix). As shown in Fig. 5, fr{d) is sam¬ 
pled in a circular annulus centered on the origin. 
The outer and inner radii of this annulus were set to 
minimize boundary effects and to achieve the asymp¬ 
totic scattering regime. The differential cross section 
da/dO —calculated by averaging \fr{d)^ across all r 
within the annulus at fixed 6 —is numerically inte¬ 
grated to produce the cross section a. 

A comparison between the total cross section a 
determined using the exact da/d6 and using du/dO 
from the Born approximation is shown in Fig. 6. The 
two approaches agree at high kinetic energy. At low 
energies, however, the Born approximation fails and 
underestimates the cross section. 

An improved scaling theory prediction for the lo¬ 
calization length ^ using the exact differential cross 
section is shown in Fig. 7. We fix n = 0.2 /im“^, 
which corresponds to a vertical slice in Fig. 4. The 
localization length predicted using the exact differ¬ 
ential cross section exceeds the Born-approximation 
prediction by more than three orders of magnitudes 
at all particle energies. Furthermore, the localiza- 



FIG. 5. Process of extracting scattering properties from 
a simulating propagation of a plane wave for the time- 
dependent Schrodinger equation. The initial wavevector 
k is in the y direction, (a) Real part of the wavefunction 
ipix, r) shown in false color. For all of the data shown in 
this figure, tk = ks x 240 nK, Vb = fcs x 1000 nK, and 
w = 400 nm. Violation of the Born approximation is 
evident as significant distortion of the initial plane wave 
by the potential barrier at the origin (not pictnred). (b) 
Real part of the scattered wavefnnction ipsdx) obtained 
by subtracting a plane wave propagated in free space 
from the data shown in (a). The color bar applies to 
the simnlated wavefnnctions shown in (a) and (b). (c) 
The scattered probability density \ipscix)\'^ is shown as 
a density plot. The cross-hatched region is sampled to 
determine that differential cross section, (d) Differen¬ 
tial cross section obtained from the data in (c). Points 
are shown at fixed 9 and different r within the circular 
annulus. 


tion length is always larger than the characteristic 
length scales of the disorder, while the Born ap¬ 
proximation predicts unphysical localization lengths 
much smaller than w and \pn « 2 /rm. For the dis¬ 
order parameters explored here and a 100 ^m system 
size, localization with be visible for point-like disor¬ 
der for particles with energies less than ks x 20 nK. 
While for a speckle potential with the same A the 
percolation threshold is /cs x 500 nK, the percola¬ 
tion threshold for this point-like disorder potential 
is approximately ks x 0.25 nK, and thus AL will be 
the dominant influence on transport. 

As a second approach to determining localiza¬ 
tion lengths, we performed complete time-dependent 
simulations of a wavepacket propagating through a 
point-like disorder potential (see Appendix B). This 
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FIG. 6. Ratio of total collision cross section determined 
from the Born approximation to cr calculated using 

the exact differential cross section for w = 400 nm and 
Vo = fcfl X 50 nK. The error bars show the spread in 
values determined across the inner and outer radii of the 
annulus displayed in Fig. 5c and include the effects of 
numerical errors and deviation from the far-held limit. 


approach is inspired by the method used to observe 
AL in ID [5, 6] and 3D [7-9] ultracold gases: atoms 
initially confined in a small region of space by a 
trap are allowed to expand into a disordered po¬ 
tential. The atoms propagate through the disor¬ 
der potential, and, through scattering, eventually 
adopt a localized profile. We numerically simulate 
independent realizations of a potential according to 
Eq. 1, with a disorder-free, circular region of radius 
i? centered on the origin (see Fig. 8a). As an ini¬ 
tial condition, a Gaussian wavepacket ip{x,t = 0) oc 
exp (ffcr — r^/2(i?/2.2)^) with kR » 1 is prepared 
in the disorder-free region. The wavenumber k is 
varied from 0.07-0.4 /rm“^, which corresponds to a 
kinetic energy from ksx 3-100 nK. The shift in the 
average kinetic energy of the wavepacket resulting 
from the gaussian envelope is less than 2% of e^. 

The wavefunction is propagated according to Eq. 9 
for total time r = 146 ms using timesteps 6t = 
14.6/iS in order to resolve the effects of the dis¬ 
order potential energy on the wavefunction. Dur¬ 
ing propagation, the wavepacket disperses into the 
disorder potential, eventually forming a localized 
state (Fig. 8b). We perform an angular average of 
\ijj{x,T)\^ at fixed radius r, which is then averaged 
over multiple realizations of the disorder potential 
to produce the disorder-averaged radial probability 
density Prir). Sample Prir) for fixed disorder pa- 



X nK) 


FIG. 7. Localization lengths predicted using scaling the¬ 
ory and the exact differential cross section (circles) and 
the Born approximation (solid line). The density for the 
disorder potential is n = 0.2 ^m“^, the average disorder 
potential energy A = fcs x 1000 nK, and the Gaussian 
barrier height Vo = ks x 9950 nK. The error bars are 
determined using the spread in the differential cross sec¬ 
tion determined across the inner and outer radii of the 
annulus displayed in Fig. 5c. 


rameters Vq = ks x 1000 nK, w = 400 nm, and 
n = 0.08 p,m“^ and kinetic energy Ck = ks x 100 nK 
are shown in Fig. 9a. An exponential distribution 
at large r is rapidly achieved, which then expands 
and relaxes at long times to a static, localized state. 
We fit Prir) to an exponential decay for r = 43- 
50 fim for ej, < /c^xlO nK and r = 57-68 /im for 
Cfe > fcsxlO nK at each r in order to determine 
the characteristic size of the propagated wave- 
function. The asymptotic localization length ^ (that 
would be observed in an experiment probing atoms 
expanding through a disordered potential) is deter¬ 
mined by fitting to an exponential function, dis¬ 
carding points at short times. 

These numerical simulations enable us to bench¬ 
mark the individual scattering technique against a 
calculation without any inherent approximations. A 
comparison is shown in Fig. 10. At the lowest ener¬ 
gies, the single scattering approximation predicts lo¬ 
calization lengths smaller than the de Broglie wave¬ 
length 27r/fc of the particle, which occurs because 
na < 27r/fc and the approximation that I a = (ncr) ^ 
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FIG. 9. Procedure to determine ^ from the time- 
dependent simulations of localization, (a) Sample radial 
probability density showing the approach to a steady- 
state profile. The propagation times are r = 3 ms 
(solid black line), 15 ms (red dashed line), 30 ms (blue 
dotted-dashed line), 90 ms (green dashed-dotted-dotted) 
line, 120 ms (short dashed line), and 150 ms (pur¬ 
ple dotted line). For the data shown in this hgure, 
tk = ks y. 100 nK, Vb = fes X lOOOnK, w = 400 nm, and 
n = 0.08 . (b) Fitted exponential decay lengths 

from data such as those in (a) for Sk = ks x 100 nK (red 
triangles), 35 nK (blue circles), and 3 nK (black squares). 
Only a hfth of the points generated for each value of tk 
are shown. The uncertainty in the points is too small to 
be visible. Fits of these data to an exponential function 
(lines) are used to extract the asymptotic value of the 
localization length 


FIG. 8. Simulated probability density distributions 
shown in false color from a time-dependent simulation of 
localization. A Gaussian wavepacket initialized at t = 0 
in the disorder-free region marked by the white dashed 
line is shown in (a). The wavefunction propagated for¬ 
ward in time for 146 ms is shown in (b). The potential 
barriers that constitute the disorder potential are mag¬ 
nified and marked in black for clarity. For these simula¬ 
tions, €k = ks y 25 nK, Vb = fcs x 1000 nK, w = 400 nm, 
and n = 0.08 


is violated. In contrast, the results of the time- 
dependent simulation always produce localization 
lengths greater than the de Broglie wavelength. As 
expected, the single scattering approximation pro¬ 
vides an upper bound on the numerically deter¬ 
mined localization length. At low energies, the sin¬ 
gle scattering approximation and results of the time- 
dependent simulations agree within an order of mag¬ 
nitude. The Born approximation, however, is four 
orders of magnitude smaller than the numerically 
simulated ^ at low energies. 

Finally, we calculate a localized density prohle for 
a thermal gas of particles under conservative con¬ 
ditions (i.e., n = 0.2/rm“^) using the simulated lo¬ 


calization lengths from Fig. 7, corresponding to the 
single scattering approximation and exact scattering 
cross section. We assume that the density profile as¬ 
sociated with each energy in the thermal ensemble is 
a radially symmetric exponential function centered 
at the origin with a decay length determined by a fit 
to the data in Fig. 7. The average density profile of 
a localized gas of 10,000 Maxwell-Boltzmann parti¬ 
cles at temperature T = 10 nK is shown in Fig. 11. 
Dimensionality plays a helpful role, since, unlike in 
three dimensions, the lowest particle energies (with 
the smallest localization lengths) are the most prob¬ 
able. Hence, most of the particles localize within a 
relatively small, experimentally accessible area. For 
example, at T = 10 nK, 90% of the particles are 
localized within a radius of 20 /rm. Furthermore, 
just 2.5% of the particles have energies below the 
percolation threshold and are classically trapped. A 
gas released into the disorder potential with a time- 
independent profile similar to that shown in Fig. 11 
is straightforward to image with high signal-to-noise 
ratio and would provide a clear signature of two- 
dimensional AL. Since the single scattering approxi¬ 
mation is an upper bound on the localization length 
observed in an expansion experiment, the localized 
thermal density profile should be smaller in extent 
than what is shown in Fig. 11. 
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Sk (nK X k^) 

FIG. 10. Comparison of localization lengths from a nu¬ 
merical simulation of the time-dependent Schrodinger 
equation (blue squares), from the perturbative correction 
to the independent scattering transport properties (red 
circles), and from the Born approximation (black line). 
For comparison, the de Broglie wavelength 2n/kis shown 
as a dotted line. The parameters used for this plot are 
Vo = ks X 1000 nK, w = 400 nm, and n = 0.08/rm~^. 
The error bars for the circles shown the impact of the 
difference in differential cross section across the inner 
and outer radii for the annulus shown in Fig. 5c. and 
the standard deviation in the energy (resulting from the 
Gaussian envelope of the wavepacket) for the squares. 



r(nm) 

FIG. 11. Thermal density profile for Vo = ks x 9947 uK, 
w = 400 nm, n = 0.2 mum~^, and T = 10 nK. Here, N is 
the two-dimensional number density for 10,000 particles. 


V. CONCLUSION 

We have identified an experimentally feasible ap¬ 
proach to observation two-dimensional AL using ul¬ 
tracold atoms using point-like disorder. Our ap¬ 
proach requires nanoKelvin-scale temperatures in 
conjunction with a dilute, strongly scattering dis¬ 
ordered potential. Using a single scattering approx¬ 
imation and exact numerical simulations, we predict 
observable localization lengths for realistic particle 
energies. This method can be extended to interact¬ 
ing systems such as a gas composed of two spin states 
of fermionic atoms, which may enable exploration of 
the analog of two-dimensional metal-insulator tran¬ 
sitions [13-15]. 

Appendix A: Individual Scattering Simulation 

In simulations of scattering from a single Gaus¬ 
sian potential barrier, square grids were employed 
ranging from — 3600 to 10000 /im^ with 
discretization sizes of 0.01 to 0.25 fim. The size of 
the time steps was scaled with the shorter of the 
two physical length scales in the system, the dis¬ 
order width w and the particle de Broglie wave¬ 
length 27r/fc, according to 6t=50.1 ns/^m. Approxi¬ 
mately 16,000 time steps were taken, which allowed 
the wavefunction initially at the center of the po¬ 
tential to scatter outward through roughly L/2, in¬ 
dependent of the group velocity. Propagation for 
longer times resulted in significant errors induced by 
boundary and finite-size effects, manifest as a peri¬ 
odic modulation of the scattered wave. The scat¬ 
tered wavefunction was sampled between radii of 3- 
5 /rm in order to determine the scattering amplitude. 

Appendix B: Time-Dependent Simulation of 
Anderson Localization 

We used a 1536 x 1536 discretization of the system 
area 349 x 349 /rm^ and populated the region r > R 
with pointlike disorder given by 

V{x)=Y,Voe-i^)\ (Bl) 

i 

where the coordinates Xi are chosen randomly with 
the constraint that Vi > R = 21.8 /im. We choose 
lu = 400 nm, Vq = ks x 1000 nK, and n = 
0.08 p.m“^, which corresponds to disorder strong 
enough to produce localization within the simulation 
space while remaining computationally tractable. 

We include an imaginary potential component 
near the edges of the system grid in order to absorb 
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probability current approaching the periodic bound¬ 
ary. This enables dynamical simulations to continue 
for the long times necessary to observe localization. 
The absorbing boundary potential is given by 


V = 



kl, \y\ > 129/rm 
otherwise 


(B2) 


with r^ax = 174 fim and i? = fcs x (5 x 10 nK. 
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